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Abstract 

A few years ago, diffraction of atoms by double slits and gratings was 
achieved for the first time, and standard optical wave-theory provided 
an excellent description of the experiments. More recently, diffraction of 
weakly bound molecules and even clusters has been observed. Due to their 
size and to possible breakup processes optical wave-theory is no longer ad- 
equate and a more sophisticated approach is needed. Moreover, surface 
effects, which can modify the diffraction pattern of atoms and molecules, 
give rise to further complications. In this article a fully quantum mechan- 
ical approach to these questions is discussed. 

1 Introduction 

New and interesting theoretical challenges are presently arising in atom and 
molecular optics, as we want to show here. The rapidly growing field of atom op- 
tics exploits wave-theoretical aspects of quantum-mechanical particles and their 
interference effects. For electrons and neutrons this has a long history, but for 
atoms it is more recent. Advances in nanotechnology have made it possible to 
manufacture structures such as double slits and gratings with slit widths as small 
as 50 nm. Typical optics experiments have been carried over to atoms. For exam- 
ple, atomic beams have been sent through a double slit and diffraction patterns 
were observed Similar experiments have been performed with transmission 
gratings and large numbers of diffraction peaks were observed with high resolu- 
tionij. 
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In such experiments the incoming atoms usually can be described by plane 
waves which are associated with quantum-mechanical point particles. Classical 
optics with its standard wave-theoretical methods and approximations, in par- 
ticular those of Huygens and Kirchhoff, has been successfully applied to atom 
optics and has yielded good agreement with experiments. 

Interesting problems with this simple picture have, however, appeared quite 
recently. One complication arises in the diffraction of molecules, in particular of 
weakly bound systems. When a diffracted system is observed at a nonzero angle 
from the incident direction, then the system has received a momentum transfer 
from the grating. For a weakly bound system this may induce breakup processes 
and result in a change of the diffraction pattern. Moreover, weakly bound systems 
are large compared to atoms, and treatment as point particles may no longer be 
adequate. In addition, low-lying molecular rotational or vibrational states may 
be excited by the passage through the grating, and this may also change the 
diffraction pattern. As a further complication, surface potentials resulting from 
van der Waals forces have been seen to change diffraction patterns [^]. 

Experimentally, molecular diffraction has been observed for He2 in Ref. and 
for Na2 in Ref. . The molecule Na2 is tightly bound and very small compared to 
present-day slit widths, and the results seem to be describable by classical wave- 
optics. The recently discovered helium dimer He2 [0, ^, however, is extremely 
weakly bound (—0.11 /leV) and quite large, with an estimated diameter of about 
6 nm . For a slit width of 50 nm this considerable size of the molecule has some 
effect on the diffraction pattern, and this effect will increase for future smaller 
slit widths. 

We have performed a theoretical analysis of diffraction of weakly bound di- 
atomic systems by a transmission grating p. The main point which we wish 
to stress in this article is that diffraction by a transmission grating really means 
scattering by the grating bars. Therefore, since breakup processes and finite-size 
effects are to be included, the proper approach is to consider quantum mechan- 
ical scattering of composite systems by an external potential, W say. For the 
beam velocities used in experiments the electronic degrees of the two atoms can 
be neglected, and thus one deals with a system of two particles, bound by a po- 
tential, denoted by V. Due to the presence of the external potential W resulting 
from the grating bars ("walls") the problem is quite different from one-particle 
scattering and rather resembles, at least mathematically, the three-body problem 
where, loosely speaking, the external potential plays the role of a third body. 

Not surprisingly one can therefore adapt methods of few-body theory to the 
present problem, in particular the Alt-Grassberger-Sandhas (AGS) equations 
which are related to the Faddeev equations. In Section 2 we briefly outline how 
one can adapt this to the problem at hand, and we obtain AGS equations for 
the scattering amplitude of a bound two-particle system by an external poten- 
tial. These equations decouple here and can be iterated in a straightforward 
way. In Section 3 we indicate how the lowest order can be specialized to diffrac- 
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tion scattering by a transmission grating with a short-range strongly repulsive 
("reflective") bar potential. Our results are illustrated for He2. At the end we 
will discuss the role of weakly attractive surface potentials and how they can be 
incorporated. 



2 Scattering of a composite two-particle system 
by an external potential 

In this section the quantum mechanical scattering problem of a bound two- 
particle system (molecule) by a quite general obstacle will be considered. The 
latter is described by an external potential of the form 

l^(Xi,X2) = W^l(Xi) + l^2(x2) (1) 

where xi 2 is the position of the respective atomic constituent and Wi^2 the re- 
spective single-atom potential. The binding potential is denoted by V{yii — X2), 
with negative binding energies and two-particle bound states (ff-yi^xi — X2), 
where the potentials V, Wi and W2 are assumed to be of short range. Among 
all possible scattering channels in the following only those will be considered in 
which the molecule is elastically scattered or changes its internal bound state. 
The plane-wave scattering matrix element for an incoming molecule with mo- 
mentum P' and internal state (py and an outgoing molecule with momentum P 



and internal state 0-^ reads |]T0 



-27ri6 ( 1^ + E^~ - t(P, 0^; P', 0y) 



(P,0,|^|P',0y) = 5,y(5(3)(P-P') (2) 

|p|2 |p/|2 
^ + E^-- 

2M ^ 2M 

where M is the molecular mass. The transition amplitudes t(P, 0^; P', 0y) de- 
termine the molecular diffraction pattern and can be calculated by adapting the 
AGS approach as follows. Introducing the Green's operators 

(3) 



Goiz) 




-^0; 




G{z) 


= 


-H^ 


-V-W) 


Gv{z) 


= 


-Ho 


-v)-\ 


Gw{z) 




-Ho 


-wr\ 



where Hq denotes the sixdimensional free Hamiltonian (kinetic energy) of two 
particles, the operators Uvv{.z) and Uwv{z) are defined by the relations 

Giz) = Gviz)+Gviz)Uvviz)Gviz), (4) 

G{z) = Gw{z)Uwv{z)Gv{z). (5) 
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As in Ref. it can be shown that Uyv is directly related to the scattering matrix 
Eq. (H) through 

t(P, (f)Y, P', 0y) = (P, (j)y\Uvv (|Pr/2M + + iO) |P', 0y). (6) 

The resolvent equation G{z) = Gv{z) + Gv{z)WG{z), applied to the left hand 
side of Eq. , gives 

Gv{z) + Gv{z)WG{z) = Gv{z) + Gv{z)Uvv{z)Gviz), (7) 

and Eq. (H) inserted into the left-hand side of Eq. (|^ gives 

Uvviz) = WGw{z)Uwv{z). (8) 

We introduce the T matrix Tw{z) = W + WGw{z)W , which satisfies 

Tw{z)Go{z) = WGw{z) {G^\z) + W) Go{z) = WGw{z). (9) 

Eq. (H) can then be written as 

Uvviz) = Twiz)Go{z)Uwv{z). (10) 

In a similar way a second AGS-like equation is obtained as 

Uwviz) = Go\z) + Tv{z)Go{z)Uvv{z), (11) 

where Tv = V + VGvV, in analogy to Tw- Once iterated Eqs. ([10|) and (0) 
decouple and give 

Uvviz) = Tw{z) + Tw{z)G,{z)Tv{z)G,{z)Uvv{z). (12) 

Although Eq. (|T^ explicitly determines only the transition operator for elastic 
and excitation processes, Uvv{z), all channels are implicitly included. Similar to 
Eq. (^) it can further be shown that 

Uoviz) = G^\z) + (1 + Tv{z)Go{z))Uvv{z) (13) 

is the transition operator for a molecular break-up into two asymptotically free 
particles. 

Eq. ([T^) is exact and contains both breakup and finite-size effects. It lends 
itself to iteration, and to lowest order in Ty one has 

Uvv = Tw (14) 

By Eq. (|]) this gives 

t(P,0^;P',0y) ^ (P,0^|T^(E + iO)|P',0y) (15) 

with E = |Pp/2M + E^. Although is the transition operator for scattering 
of two asymptotically free particles by the external potential W it is not on the 
energy shell in Eq. (p^. The interaction V is taken into account in Eq. (pTSp 
through the wave-functions 0^ and (py. 
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3 Application to diffraction gratings 

Although the potential W = Wi + W2 is additive in the two particles and looks 
very simple, Eq. ( ]T^ ) is not straightforward to evaluate since the Green's operator 
Gw is not easily expressed by one-particle Green's operators. In the diffraction 
domain the de Broglie wavelength is small compared to the extension of the 
obstacle and the potential. To evaluate Eq. (|T5[) for small scattering angles one 
obtains an integral expansion of Eq. (plSf ) by inserting two complete sets of two- 
particle plane waves, |pi,P2) and |p']^,p2), on the right-hand-side. The kinetic 
energies of these states are not equal to E (i.e. Tw is not on the energy shell). 
However, the momentum space wave- functions 0^(p) and (p) are centered at 
relative momenta p which are small compared to P = |P|, and this fact can be 
used to show that in Tiy(i? + iO) one can replace E by pf/2m + p2/2m, where, for 
simplicity, the atomic constituents are assumed to have the same mass, m = M/2. 
The Lippmann-Schwinger equation for the potential W further enables one to 
reduce Eq. (|T5|) to one-particle transition amplitudes for the potentials Wi and 
W2- For reflecting bars these can be calculated by standard methods. 
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Figure 1: Transmission grating with slit width s and grating period d. 

If the bar and slit widths are much larger than the diameter of the two-particle 
system, then the one-particle amplitudes are sharply peaked at the forward direc- 
tion on a scale on which the other expressions appearing in the integral expansion 
of Eq. (|15D are slowly varying and can therefore be dropped. This then leads to 
the usual diffraction pattern of a point particle of mass 2m. 
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When bar and slit widths decrease the situation changes. In this case one 
obtains a more comphcated coherent amphtude part, tcoh, the superposition of 
amphtudes for the interaction of the molecule with individual bars. The incoher- 
ent part from the interaction with several bars becomes appreciable only when d 
and s (see Fig. ^ are of the order of the molecular diameter. The net result is of 
the form 

^C'(osP/M;P.)5«f.(P,) (16) 

where the first factor is the molecular transition amplitude for a single bar of 
width d — s, for total momentum P and total mass M. The second factor is 
the usual sharply peaked grating function known from optics, and it gives the 
diffraction peaks at the same locations as for point-particles with corresponding 
lateral momentum transfer. The S function is due to the assumed infinite extent 
of the grating in the third direction. The amplitude of the diffraction peaks 
is determined by t^°^. The latter has been explicitly calculated by means of 
Eq. (|T3]) in Ref. |^ and it depends only on the absolute value of the ground-state 
wave-function. One obtains 

Cf(7,P/M;P2) = - (^^M -'"^^-^^-^/^^^ / d% 2 e^^^-^/^^|0,(x) P (17) 
+ / dxidxa J^~' dx2|0T,(x)|2 sin [^{d - s - xa)] /P2 ■ 

The expression preceding the first integral is the single-bar amplitude for a point 
particle of mass M and total momentum P. The complete expression on the r.h.s. 
reduces to this if |(/)-),(x)p contracts to a point. 

In Fig. I we have plotted P for He2 and the corresponding quantity for 
a point particle. The bar width is 25 nm. For a symmetric grating of period 
50 nm the vertical lines in the inset indicate where the peaks of the grating 
functions cut out the diffraction peaks of first order, second order and so on. For 
a symmetric grating and point particles there are no even-order diffraction peaks 
since the single-bar amplitude vanishes there, as seen in the inset. For the dimer 
this is not true, and it therefore has small even-order peaks which become more 
pronounced for smaller bar and slit widths. 

Intuitively we expect that finite-size effects of He2 might be partially taken into 
account by a larger bar and smaller slit width. Fig. |^ shows that the single-bar 
amplitude of He2, for a bar width of 25 nm, can be approximated for small lateral 
momentum transfer by that of a point particle for bar width of (25-1-2.8) nm. For 
a grating this would mean a nonsymmetric grating of correspondingly smaller slit 
width. The expectation value of \x2\ for the ground-state wave-function of He2 
we are using is just 2.8 nm. 
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Figure 2: Diffraction by a single bar of width 25 nm for lielium dimers and point 
particles. For a symmetric grating with this bar and slit width the vertical lines 
in the inset pick out the corresponding diffraction peaks. For dimers there are 
even order peaks which are absent for point particles. The lateral momentum 
transfer is P2 = T1K2. 

4 The role of attractive surface potentials 

When a rare-gas atom such as helium passes through a slit it experiences a 
repulsive potential of very short range. However, if the atom is polarizable there 
will also be a longer-range attractive surface potential acting on it, due to van der 
Waals forces. Recent experiments have shown that this additional attractive 
potential can have an unexpectedly large influence on the diffraction pattern, not 
so much so for helium but rather for the heavier rare-gas atoms, which have a 
higher polarizability. 

The attractive part of the surface potential, denoted by Ws-, is of the form 



where is the distance of the atom from the slit surface and where C3 is a 
constant depending on the grating material and the particular rare-gas atom. 
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Figure 3: Diffraction of helium dimers (solid line) by a single bar of width 25 nm 
and diffraction of point particles (dashed line) by a single bar of width (25 + 2.8) 
nm. There is agreement for small lateral momentum transfer P2 = T1K2. 

In a first step we have taken Ws into account for the atomic case, with quite 
an unexpected and surprising result. The simple picture of the grating given in 
Fig. 1^ as an array of rectangular bars is not adequate if the surface potential 
plays a role. In fact, with the simple geometry of Fig. ^ our results did not even 
qualitatively agree with early experiments for heavier rare-gas atoms. As it then 
turned out, the etching process produces bars whose cross section is not quite 
rectangular but rather trapeze-like as in Fig. |, with a small angle a of about 
8° off the perpendicular direction. This seems to be a small deviation from the 
idealized rectangular form of Fig. |l|, but taking the angle a into account gave at 
least qualitative agreement. In the mean time the approach, both experimentally 
and theoretically, has been so perfected, that it has been possible to determine 
both the grating parameters as well as the surface potential constant C3 from 
experimental data 0]. 

For dimers the surface potential can be taken into account by adding W^^ 
and to W in Eq. (]^). The details are presently being worked out. 
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Figure 4: Trapeze-like cross section of grating bars. 

5 Conclusions 

We have shown in this article that new and interesting theoretical questions arise 
in atom optics when the spatial extent of a diffracted system can no longer be 
ignored and when possible excitation energies become comparable to the kinetic 
energy of diffracted systems. Here we have studied the diffraction of a weakly 
bound two-particle system by an obstacle, in particular by a transmission grat- 
ing. We have obtained an analytic expression for the elastic diffraction amplitude 
which incorporates breakup processes and finite-size effects. This expression con- 
tains the familiar grating function of classical optics as a factor and this puts the 
diffraction peaks at the angles expected from classical wave theory. This grating 
function is multiplied by the molecular diffraction amplitude of a single grating 
bar, and in this factor a difference from the point-particle case may arise. 

For a molecular diameter very small compared to the bar and slit widths 
the derivations from the point-particle result are indeed negligible, as expected. 
This is no longer true for smaller bar and slit widths. In particular, heights 
of diffraction peaks may become lower than for point-particles. This can be 
attributed to breakups that diminish the actual number of molecules in a given 
direction. Moreover, the decrease in peak height goes with the diffraction angle 
and thus with the lateral momentum transfer, which is quite reasonable since 
higher momentum transfer should imply more breakups. In addition, diffraction 
orders that are suppressed to zero for point particles and for symmetric gratings 
may reappear. 
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If the slits become smaller than the diameter of the two-particle system, ad- 
ditional terms appear in the diffraction patterns. But this case is presently not 
of experimental interest. 

The appearance of the diffraction peaks at the wave-theoretical locations for 
gratings holds quite generally, also for larger composite systems. For dimers a 



simplified argument has been given by us in Ref. |rT[, and the general case will 
be treated elsewhere. 
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